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We determine numerically the ground state of the two-dimensional, fully polarized electron gas 
within the Hartree-Fock approximation without imposing any particular symmetries on the solu- 
tions. At low electronic densities, the Wigner crystal solution is stable, but for higher densities (r s 
less than ~ 2.7) we obtain a ground state of different symmetry: the charge density forms a trian- 
gular lattice with about 11% more sites than electrons. We prove analytically that this conducting 
state with broken translational symmetry has lower energy than the uniform Fermi gas state in the 
high density region giving rise to a metal to insulator transition. 

I. INTRODUCTION 

The two-dimensional homogeneous electron gas is one of the fundamental models in condensed matter physics. 
Despite its simplicity - the system consists of electrons interacting through a 1/r-potential to which a unif orm 
positive background is added for charge neutrality - the phase diagr am at zero temperature is nontrivia l 1 | 2 | 3 i In 
general, it is given in terms of the dimensionless parameter r s = 1/ ^/nna 2 B , where n is the electronic density and 
olb = h 2 /(me 2 ) the Bohr radius (see Section [iv] for notations and units). At low density (large r s ), the potential 
energy dominates over the kinetic energy and the system forms a triangular lattice, the Wigner crystal (WC), whereas 
in the high density region (r s — > 0) the kinetic energy favors a uniform Fermi gas (FG) phase^ (the simplest state 
given as the determinant of the plane waves with wave vector k of modulus smaller than kp). The energy of the 
FG is known analytically. Already Wigner 4 argued that the unpolarized FG is unstable even in the limit r s — » 0. 
Later, Overhauser claimed the instability of the unpolarized WC with respect to spin-density waves, even within 
the Hartree-Fock approximation (HF) 5 . It has further been conjectured that the Coulomb potential prevents any 
first order transition between the WC and a FG 6 . Despite these rather general instability theorems, there are few 
quantitative calculations of the true ground state of the electron gas within HF^. A previous HF study 8 of the two 
and three dimensional electron gas compares the FG energy with the energy of various states with imposed crystal 
symmetries. For the polarized two-dimensional gas, they find lower energies for a crystal for r s larger than 2. Only 
recently an unrestricted HF study of the unpolarized three-dimensional electron gas was performed which proposes a 
more complicated structure of a ground state with spin-density waves in the high density region 9 . 

Indeed, establishing the precise HF phase diagram of the electron gas influences the correlation energy estimations, 
since by definition the many-body correlation effects must be evaluated with respect to the true HF ground state. 
Further, even in more advanced technics, the antisymmetry of the wavefunctions is in general provided by a single 
Slater determinant. 

In the present study, we consider the fully polarized two-dimensional electron gas, analytically, and numerically. 
Section O summarized our numerical results. At low densities, r s > 2.7, our simulations always lead to a WC. 
For higher densities and large enough number of electrons, TV, the solution is neither a FG nor a WC: the density 
modulation corresponds to a partially occupied crystal of different symmetry compared to the WC phase. As the 
number of sites is larger than A/", we refer to this solution as a metallic phase. Details of our numerical methods are 
given in sect ion |III| 

In Section |IV| we remind some definitions and notations particularly used in the following. Section [V] is devoted to 
derive rigorous, analytical upper bounds on the energy of the metallic phase. These bounds are obtained in the limit 
r s — > where the calculation is simplified by the long range behavior of the interaction potential. In the conclusion, 
Section |Vl} we briefly summarize the results of the paper and discuss their relevance. 

The new HF solutions discussed in this paper within the HF approximation open a new perspective for the qualitative 
understanding of the experimental observed metal to insulator transition^ and should be considered in studies beyond 
the HF approximation. The possibility that the experimental findings are driven by interaction effects - and not by 
disorder - was recently adressed in Ref J 11 * 12 1 considering an extended Hubbard model. 
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FIG. 1: Extrapolated energies Eoo — Efg (milli-Hartree units) versus r s . Points with error bars: present calculations, full line 
(red): data of RefP, full line (green): fit to present results for (r a < 2.7), vertical dash line (blue): r s ~ 2.7. 



II. NUMERICAL RESULTS 

The N-body Hamiltonian, H = K + V, contains the kinetic energy K and the 1/r-periodic Coulomb potential 
V where a uniform positive charge background is subtracted. Within the HF approximation, the search of the true 
ground state of the quantum many-body system is reduced to the simpler problem of finding the lowest energy states 
in the subset of the Slater determinants (see Eq. |9|. Let <E> = A • • • A ^at be the Slater determinant associated with 
the single particle states {^i} an d "02? ■ ■ ■ ? ^n) the corresponding energy expectation value. 

By a kind of descent method described in the following section, we numerically study systems with up to 500 electrons 
at densities corresponding to r s = 1 up to r s = 30. Since we expect the electrons to crystallize on a triangular lattice 
at low densities, we choose periodic conditions compatible with this geometry. Thus, the unit cell of the periodized 
system is given in terms of two vectors {Li, L2} of length L and with an angle of 60 degrees between both; the volume 
of the unit cell is Vt = L 2 ^/3/2. We have restricted our study to system sizes which are compatible both with the 
triangular lattice and with a closed shell occupation in /c-space. Any triangular crystal with unit-cell vectors {ei,e2} 
is compatible with the boundary conditions if it satisfies Li = lei -\-rne2 and L2 = — raei + (Z + ra)e2, where (Z, ra) are 
two non-negative integers. The number of sites of the lattice is given by N c = det(Li, L 2 )/ det(ei, e 2 ) = Z 2 -f ra 2 -f Zra. 

In Figure [l] we report the energies of the obtained HF ground state £"00 (r s ), extrapolated to the thermodynamic 
limit, as a function of r s . In the low density region, for r s > 3, we obtain good agreement with the results of Trail et 
aZP, which imposed a ground state build as a complete band of Bloch wave functions of the triangular WC lattice. 
But for smaller r s we find lower energies which remain also below the FG energy down to r s = 1. 

The Figure [2] (a) shows the typical charge density for r s ~ 3. In this case we have a triangular lattice with exactly 
N sites. The Figure [2] (c) shows the Fourier transform of the charge density. The support of the Fourier transform is 
the six-fold star corresponding to the triangular lattice of the charge density with N c = N. 

For 1 < r s < 2.7, we find new kind of ground states (see Figure [2] (b)): the support of the Fourier transform ((see 
Figure [2] (d)) is still a six- fold star corresponding to a triangular lattice, but this triangular lattice has a number 
of sites N c larger than N. The system lowers its energy by delocalizing the electrons on a denser lattice with more 
sites than electrons (Figure [2]). This denser lattice is characterized by integral numbers (Z / ,ra/) different from the WC 
lattice (Z, ra). For some system sizes TV, the maxima of the Fourier transform correspond to various couples of {Z', ra/} 
leading to different number of lattice sites, N c = V 2 + ra/ 2 + I'm! . In any case, the system looks like a periodic crystal 
with an incomplete band in contrast to the WC solutions of fully occupied bands, studied in Ref. 8 . We refer to this 
solution as the metallic phase. However, as r s approaches zero, the energy gain of this metallic crystal compared to 
the FG gets more and more tiny. At the same time, N c is either constant or increases when r s decreases (apart for a 
few exceptions). At r s < 1, the FG solution is stable for our finite system sizes (N < 500). 

Now we would like to understand the nature of the Slater determinants in the metallic phase. A Slater determinant 
is obtained as a set of N orthonormal single particle wavefunctions Only the space generated by the ^'s is 
relevant, and in order to understand the numerical results we need to choose a canonical representation of the ^'s. 

Let {4>i}i=i...N be a basis corresponding to some indexation of the plane waves associated to the wave vectors of 
the Fermi sphere. As r s is small, the space generated by the ^'s becomes close to the space generated by the 0's. Let 
M be the square matrix defined by = {(j>i\^j) which measures the overlap of the two Slater determinants. The 
Singular Value Decomposition (SVD) of M is M = UcrV where U and V are unitary matrices and a is a diagonal 
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FIG. 2: Left : charge density p(r)/(p) — 1, with p(r) = J^iLi l^( r )| 2 - Right: p(k), the Fourier transform of the charge density, 
where p(0) = iV has been removed. The grid-points are thoose compatible with the periodic conditions. Top: the number of 
maxima is N c = 9 2 + 5 2 + 9 x 5 = 151. Bottom: the number of maxima is N c = 12 2 + l 2 + 12 x 1 = 157. Gray levels corresponds 
to the same density in both figures. Colored lines correspond to Li = lei +me2, where the numbers stands for /, m (see text). 



positive matrix. Then the orthonormal set {^}i=i...iv defined by: 

^E^fc (!) 

k 

is a basis of Span({^}) close to the basis {<f>i}i=i...N of Span({^}) in the sense that ((^|?//) = (UaU T )ij « 5^ as 
soon as a^s is close to one. 

From now on, we assume that the single particle wavefunctions, have been chosen in this way. 

Thus if r s is not too large, ipi is close to fa (at least for i associated to a wave vector not too close to the Fermi 
surface) that is ^(k^) is close to one. Thus, the largest amplitude of ^(k) is for k = k^, and Fig J3] represents the 
next largest amplitude of that we denote &k i? fc> r 499 electrons at r s = 2.7 in 2D. 

For |k^| « hp, fa has essentially only two nonzero components: one at k = k^ and the other one at the vector 
k = k^ + where is the vector of the six- fold star of Figure 2] (d) such that k is close to the Fermi surface. This 
condition can only be satisfied for a set of k^ closed to a six-fold star as we see on Figure [3J 

One can understand why metallic states should exist at small r s in the thermodynamic limit. Let us replace a 
plane wave state k of the FG (||k|| < kp) by a superposition of two plane waves with wavevectors k and k + q 
(||k + q|| > kp). Choosing q on the six- fold star of a triangular lattice we certainly obtain a gain in potential energy. 
The increase of kinetic energy is minimized if ||k|| ~ kp and ||k + q|| ~ hp. Then, the number of solutions for k 
is optimal if ||q|| ~ 2k p. This solution corresponds to a triangular lattice of length L c = 2tt / f (\^3kp) in real space 
leading to a unit cell of volume ft c = Since the system is contained in the volume ft = \/ZL 2 /2, we will 

obtain N c = Q c /Q lattice sites, or 

N c = —N « 1.17V, (2) 

7T 

where we have used n = N/ ft = k 2 F / (Air) . 

We will elaborate this argument into an analytical proof of an upper bound on the ground state energy in section 

El 
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FIG. 3: Numerical values of bk for 499 electrons in two dimensions at r s — 2. 



III. THE DESCENT METHOD 

In this section we provide some details about the descent method used to obtain numerically the HF ground states 
of the electron gas. The variation of the total energy £(ipi,ip2, • • • ? ^n) with respect to a variation of the single particle 
state Sipi is given by 

se = ^(/i*Vi|^i> + ^Wil^)- (3) 

i i 

where hy, the so-called HF Hamiltonian, is a single particle operator depending on the full state ^ (not on the 
particular choice of the ipis). Extremal states must satisfy the following equation 

hyifri = y^Cijipj, (4) 

3 

where Cij are the Lagrange coefficients associated with the normalization constraint = Sij- Conversely, if 

f = A ■• • A ipN is not an extremum, we have 

3 

where the 0^s satisfy (^i\0j) = 0, Vi, j. Within the steepest descent method one chooses first a TV x N unitary 
transformation A = (a^) such that one obtains 

<$|0j) = O, oc <^-> = Sij, Vi,j (6) 

for the transformed single particle states 0[ — . aijOj, = dijipj. The energy £(^i + A0i, ?Pn + A#tv) can be 
expressed as a sum of rational fractions whose numerators and denominators are polynomials of order four, at most. 
Thus, it is possible to find the best A and to iterate the process until a stationary state is reached. 

In fact, this method has the same drawbacks as the steepest decent method in linear optimization problems; in 
general, it converges slowly. For linear problems, conjugate gradient methods are preferabl e 13 * 14 ^. However, since the 
HF states do not form a linear space, the genuine conjugate gradient method does not apply here. We have therefore 
adapted a variant of this method to the non-linear case. Let rji be the previous variation Sifti, and 6i is obtained 
by Eq. We then compute + A#i + \±r\\, • • • , ip n + XOn + I^Vn) fc> r six values of the pairs {A, /i} in order to 
approximate £ by a polynomial of order two in A and \i. Minimizing the polynomial with respect to A and /i, we obtain 
the new changes of the single particle states, 5ipi, and the corresponding energy change. This process is iterated until 
the relative variation of the energy, 5£/£, is sufficiently small. 

We compute the wavefunction on a N g x N g grid, the fast Fourier transform is used to switch between real and 
reciprocal spaced We have systematically checked the convergence of the solution with respect to the grid size. For 
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the FG ground state, convergence is reached once all k- vectors up to 2k p are represented in the grid (N g ~ A^N/tt). 
At larger r s , in the WC phase, the wave functions are essentially Gaussian^. The width 5 of the Gaussians scales 
as S/L oc (r s TV) -1 / 2 . For a correct resolution of the Gaussians we need L/N g oc 5, so that the number of grid points 
increases at low densities, N g oc (7Vr s ) 1//2 . Convergence is reached for N g = 32 (resp. 64, 128) for N < 43 (resp. 
N < 200, N < 500) up to r s = 30. Whenever the number of grid points is chosen too small, solutions without any 
particular symmetries are obtained. 

We have further studied the influence of the initial state on the final solution, by choosing different types of 
wavefunction for initialization: a WC state, a converged state stored at larger or lower r s , a state initialized with 
random numbers, or a "metallic state" as described above. 

Typically, the energies decrease exponentially with the number of iterations. The decrease in energy during transi- 
tions to a different symmetry is in general much smaller than the convergence within the same symmetry. We have 
often seen energy plateaus with changes of relative energy < 10 -4 just before the occurrence of a transition to a 
completely different state. For system sizes up to N = 151, the minimization is continued until a relative precision of 
10 -12 is reached, and for larger TV a relative precision of 10 -5 is used. 



IV. ENERGY OF THE POLARIZED ELECTRON GAS FOR A SLATER STATE 

In this section we set our notations and recall the basic formulas of the electron gas. We consider the Hamiltonian 
of N electrons in a 2D or 3D square box of volume Q with periodic boundary conditions: 

h 2 e 2 

H = - — A+-V (7) 
2m 2 v J 

where V is the 2-body Coulomb potential Vl r « — r j I? the electron mass is m, and e is its charge. It is convenient 

to choose Hartree as the unit of energy, Ha = ft 2 /(ma 2 B ). We get: 

H= a f(-A + ^V) (8) 

Let ip n be an orthonormalized set of N vectors of L 2 (Q). They define the TV-particle Slater determinant \I/ = f\ n ip n . 
And the energy of ^ is: 

£={*\H\9) = ?f (-^<VgA|Vg + ^^<^AVvWnAVv> J (9) 

\ n n,n' I 

where v is defined as: 

(<pi <8><p2|*#i ^^2) = / dxdyTpi{x)Tp^{y)- — 1 — rr^i 0*0^2 (z/)- (10) 

J If -2/11 

In order to avoid problems due to the Coulomb singularity, we introduce the jellium model and define the potential 
acting on the plane waves as: 

7T / 2 \ D ~ 1 

{<t>k ® <t>W \v\<j>k+q ® <j>k>-q) = ^ ^ ) (11) 

for q 7^ and otherwise, so that the total charge of the electrons is compensated by a positive background charge. 
The Fermi gas is defined by $ = A|fc|<fc F ^k where (aokp) — (27r) D N/Q and a£ is the volume of the unit sphere. 

\\k\<k F a |fc|,|fc'|<fc F 1 1 / 

As Q goes to 00 with VL/N fixed, the thermodynamic limit for the energy per particle is obtained by the substitution 



£ FG _a% n ( f 2 1 f i_ 



Q ; k°+ 2 If dkk 2 - - - * / dkdk' - l m \ (13) 
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FIG. 4: The circle is the Fermi surface. The shaded surfaces are the regions where b(k) is nonzero. The new state ipk mixing 
(j)k and 0fc+Q fc is now resonant with ip k t. 



From the definition of r s — (o^a^n 1 ^) 1 and kp-, it follows that kpa^r s clb — 27r. Thus, we have: 

2?f2 dk k 2 - T -f£ [ dkdk' j-rjr-T ) (14) 



N ag+Vf \J W<1 ** D J\ k \, W \<i ' \k -k'\D-i 

which gives for D = 2 (a| = 



If dkk 2 -^ f dkdk' — ] (15) 



V. HARTREE-FOCK UPPER BOUNDS FOR THE POLARIZED 2D ELECTRON GAS 

In this section we estimate the energy for a class of states inspired by our numerical results. 
Let us consider a state ^ = A|fc|<fc F ^fe where: 

ij)k = a k (\) k + bk(l)k+Q k (16) 

with Qk in {— 2&i?(cosp7r/3, sinp7r/3)} p =o...5- For /c = |fc|(cos sin#) we choose Qk such that |fc + Qfe| is minimal; 
that is, we choose p as the integer part of (30 /tt + 1/2) and we must assume bk is zero if k is zero or 6 = tt/6 + n7r/3. 

Furthermore, we assume that a& and bk are real positive number and invariant thru the rotation of 2n7r/6 and 
the symmetry 0^—0 (i.e., the dihedral group Dq). The t/^'s are normalized, so that a 2 , + b\ — 1 and 6^ = if 
\k - Qk\ < 2k F (l — e) (i.e., bk is not zero only in the vicinity of {kp(cosp7r/3, sinp7r/3)} p =o...5), see FigjZJ 

Thus, from Eqs |9| and [TH], the limit energy per particle is given by: 

4 = A ( / dfe (^|-A|^) + ^ / dkdk' ^-(^A^H^A^)) (17) 



where, as in (15), the fc's have been renormalized by kp and thus \Qk\ — 2. 
We define AE by: 



S ~ £fG \AE (18) 



Then 



A£ = / dk [m - A|^ fc ) - fe 2 ] + ^A£ v 



(19) 
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where 



AE V = / dkdk' (fa A fa>\v\fa A Vv> + l7 1 y/ , ) (20) 



A. Potential energy contribution: AEy 

Setting v q = l/|g|: 

^ (^fe A ^/c'M^fc A t/v) + v k -k' = (vfe-fc/ - v k - k >- Qh ,)b\,a\ 



Eq.(|20|), may be divided into 4 parts: 

• = 0, 6/c = 0}: the contribution is zero. 

• {bk' = 0, b k 7^ 0}, {6fe/ 7^ 0, 6fc = 0}: both cases are equivalent. 

is: 



+ (v k -k' - v k+Qk - k ,)b 2 k a 2 k , 

+ K-fe' - v k+Qk -Q kl - k >)blb 2 k/ 

— 2v k - k 'a k b k a k 'b k '5Q k -Q k , 

~ ( v k+Q k -k' + Vk-k'-Q k , )aka k >b k b k >SQ k +Q k/ (21) 



For = 0, 6 fe 7^ 0}, the integrant of Eq. 20 



^feA^,^ ^Afe)+ | fc j; fc/ | = K-fe' -Vk+Q k -k') b l ( 22 ) 

Let So be the sector of unit disk between — 7r/6 and 7r/6 (see Fig. |4|; then in this sector Q k = (—2,0) and by 
symmetry: 



/ dkdk' (v k -k f ~ v k+Q k -k>)b 2 k = 6 / dkdk' (v k -k f ~ Vk+Q k -k>)b 

Jb k ,=0 JkeS o ,b k ,=0 



(23) 



= 6 / dkdk' (v k -k> ~ v~ k _ k ,)b 2 k (24) 

JkeS o ,b k ,=0 

< Ce 3 + 6 / dfedfe 7 - Vjfe_ fe /)6fc (25) 

^/eG5 ,|/c^|<l-e 

where fc = (2 — k X: k y ). In So, fc = (k X: k y ) where k x is close to 1 and setting k x = 1 — x : we assume from now 
that = b(x/e). 

In Appendix A, we prove that: 

/ dfcdfc' (v k - k > - v k+Qk _ k ,)b 2 k < 8e 2 V2~e [lne" 1 + 0(1)] f dxb 2 (x)x^x (26) 

JkeS ,\k' x \<l-e JO 

{h> ^ 0, fr fc ^ 0}: 

By symmetry we can assume that k belongs to So- If k' So U S3 all the u appearing in (21) are uniformly 
bounded. And since the k-volume for each sector goes like ey/e, the contribution of these terms is bounded by 
Ce 3 . In the same way v k - k > is bounded when k' G S3 and v k +Q k - k > is bounded when k' G So- Thus setting: 

/ := ajfebjfe, + b 2 k a\, - 2a k a k >b k b k > = (a k b k > - b k a k >) 2 (27) 

g := a 2 k b 2 k , + 6^a^, + 2a k a k >b k b k > = (a fe 6fe' + &fcafc') 2 ( 28 ) 

one can check that: 

/ i~ W>fc A^fc'l^l^/c A +v k -k f < Ce 3 + 6 / dfedfe' (v k - k > f - v k+Qk+k >g) (29) 

Jb k ,b k ,^o Z7T Jk,k'es 



In Appendix B we prove that 



/ dkdkf (vk-k f f ~ Vk+Q k +k>g) < 4e 2 v / 2e [lne 1 + 0(1)] / dxyfx I dx' (f(ex, ex') - g(ex, ex' 

Jk,k'eS Jo Jx 



)) (30) 



Thus, summing the four contribution gives: 



AE V < Ce 3 + 6e 2 V2e~ [lne -1 +0(1)] dx^fx (^16b 2 (x)x + 4 J dx (/(a?, x f ) - g(x, x'))j (31) 



B. Kinetic energy contribution: 

The variation of the kinetic energy is given by: 

dk 



[ dk [(^ k \ - A\il> k ) - k 2 } = 6 / dk((^ k \ - A\^j k ) - k 2 ) 

J\k\<l JkeS 



= 6 / dx 2y m 4xb (x/e) 
Jo 

< 6 x 8e 2 V2e dx ^x xb 2 (x) 
Jo 

C. Total energy: 

Inserting Eqs ( 31|34 ) in Eq. (19), the variation of the total energy from the Fermi gas energy becomes: 
AE < 6e 2 V2e dx^fx (%xb 2 {x) + ^ [in e" 1 + 0(1)] ^166 2 (x)x + 4 J dx' (f(x, x f ) - g(x, x f )) 
= 6x8e 2 v / 2e^ dx^x (^xb 2 (x) + ^ [lne -1 + 0(1)] (b 2 (x)x - a(x)b(x) J dx f a{x f )b{x , )) S j 



Let us set 



S = e 2 V~e 



h 



I dx^fxxb 2 {x) 
Jo 



5tt 



dx\fx ( — b 2 (x)x + a(x)b(x) / dx'a(x')b{x' 



Then 



AE < 6 x SV2S [h - r s I 2 (\n S' 1 + 0(1))] 
If I2 > 0, as r s goes to 0, AE is minimal in Eq. (l39|) for S defined by: 



r 1 ( h 

<3min = -exp 

e V hr.c 



and finally inserting 8 m i n in Eq. (39) gives: 



AJ?< 6x8^/2 / 7i 

A£ < exp -- — r s I 2 

e V 



(32) 
(33) 
(34) 



(35) 

(36) 
(37) 

(38) 
(39) 
(40) 

(41) 



We now have to find a solution b(x) such that I2 is positive. Choosing b{x) = bo or b(x) =b$(l — x) leads to negative 
12. In the Appendix C, as r s goes to we find a family of b leading to : 



AE < -r s exp 



5tt 0(1) 



3r« 



(42) 
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FIG. 5: Effect of the screening on the ground state phase diagram for N = 61, 151, and 499. The different phases are labelled 
"FG" ( Fermi-gas ground state), "M" (metallic ground state), and "W" ( Wigner crystal). For all values of the screening 
parameter a of the screened Coulomb potential V(q), we have found a metallic phase of energy lower than the Fermi gas and 
the Wigner crystal energy. 



Though such a bound is correct in the thermodynamic limit for r s — > 0, this behavior in not so relevant for finite 
systems. Our numerical calculations of section [□] consider about 10 3 electrons where the uniform Fermi gas remains 



the ground state for r s < 1. Thus the asymptotic bound (|42| is not very helpful in comparing with our numerical 
results obtained for r s « 1. 

Nevertheless, for finite r s , on can choose a suitable function b and evaluate numerically 1\ and I2. For instance, 
with b = brj as in (72 ) of Appendix C and 77 = 0.001 we get 



£ — £fg 
N 



< -2. x HTV^exp 



18.5 



(43) 



We can also understand why the metallic phase does not neccessarily appear in small sized systems. For finite 
systems of N electrons, one must have at least one plane wave in the shaded region of figure |4j \k • Qk\ > 2k%(l — e). 
This gives the condition Ne^/e > 1 and using Eqs ( 36|40 ), this leads N > exp(3/i/5/2^s)- Analogous to Eq. (42), we 
obtain the following lower bound 



N > cxp ( ^ 



(44) 



i.e. N > 500 for r s = 1.8. 
disappears at r s = 1 for N 
numerical calculations. 



This bound is compatible with our numerical simulations where the metallic phase 
= 500, and may explain why the metallic phase has not been observed in previous 



VI. CONCLUSION 



Using a descent algorithm, we have computed the ground state of up to N = 500 electrons. For 1 < r s < 3, our 
solutions have lower energies than the FG or WC. These solutions correspond to denser lattices that the WC solutions, 
that is with less than one electron per site as in a metallic material. 

We have proven, in the thermodynamic limit, that for sufficiently small r s , these metallic states have always a 
smaller energy than the Fermi gas. To our knowledge, it is the first time that rigorous upper bounds for the ground 
state energy of the polarized electron gas are obtained, demonstrating that the FG is not the ground state even at 
small r s . 

Our proof relies on the behavior at infinity of the Coulomb potential, so it may be interesting to check the existence 
of these states in the case of a screened Coulomb potential. A rigorous extension of our proof is not straightforward. 
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However, as shown in Fig. [5j our numerical calculations indicate that the metallic phase persists in the presence 
of screening, at least for various system sizes and screening parameters studied. Thus, such metallic states should 
be considered as relevant candidates for further studies beyond the Hartree-Fock approximation, since, qualitatively, 
correlation effects amounts to an effective screening of the electron interaction in the high density limit. 

VII. APPENDIX A 

We have to estimate 

/(/) = ji*n,iifc'n<i dkdk ' {vk - k ' ~ "s-k'M 1 - k ^ ^ 

l-k x <e,\k' x \<l-e 

where k = (2 — k x , k y ) and / is a positive function. 

/ dk' (v k . k , - v- k _ k> ) = f dk' x asinh f^f + asinh ^±h. _ asinh V'm-ky _ ^ v' m + k (4g) 



where y r m — J I — k' x . And since asinh x — asinh y < \nx/y for x > y > 0: 



dk' (v k . k , - v k _ k> ) < I 1 £ dk' x 2 In 2 , kx _ ,, K (47) 



i+e k x k' x 



We set k x = 1 — x and y rn = \Jlx — x 2 , 



1(f) < I dxf(x) 2y m [ dk' x 2\n\ 

JO J-l+e 1 

= 4 / dxf(x) y m \ du In 

J0 Je U ~ 



x - k' 



X 



X 

2-e 



(48) 
(49) 



< 4 / dxf(x) y m I du (50) 

JO Je U — X 

= 8e 2 V^ [lne -1 + 0(1)] / dxf(ex)x^x (51) 

Jo 



VIII. APPENDIX B 



We have to estimate 



J(/,5) = J\\kUk'\\<i dMk ' ( Vk ~ k 'f - V 'k+k'9) 



(52) 



l-k x A-k'<e 



52 



can be rewritten: 



where k = (k x — 2, k y ) and / and g are positive functions of 1 — k x and 1 — k x > . Setting k x = 1 — k' x = 1 — 
k y =y, k 'y = V', and r ± = \/{x±x') 2 + (y-y') 2 , Eq. 

i(f,g)-- 



dx I dx' I dydy' ( — / < 

2 / dx I dx' I dydy ( — / < 



(53) 
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where y and y' must satisfy (1 — x) 2 + y 2 < 1 and (1 — x') 2 + ?/ 2 < 1. 
Since asinh a; < ln2(x + 1), the first term in 53 is bounded by: 

2 [ dx [ da/ f dydy' —f = 2 f dx [ dx'f F™ dy (asinh y ™ + V + asinh (54) 

JO Jx J r - Jo Jx J-Vm X — X X — X 

<2 f dx [ dx'f f Vrn dy 2 asinh (55) 

Jo J x J —y m X — X 



X — X 



;) (56) 



<4 / da; / dx'f 2y m ln(2 

JO Jx 

<4e 2 v / 2^[ln(e- 1 ) + 0(l)] / dz / dx f(ex,ex')y/x (57) 

JO Ja? 



On the other hand, using asinh a: > ln2x, the last term of (53) is: 

2 [ dx [ dx'g [ dydy'— = 2 [ dx f dx' g F™ dy asinh ^ll + asinh ( 58 ) 

JO Jx J r + JO Jx J -ym X -\- X X + X 



>4 / da? / da/#2/ m [in e" 1 +0(1)] 

JO Jx 



I(f,g) < Ae 2 V2~e [lne" 1 + 0(1)] f dx^x [ dx' (f(ex, ex') - g(ex, ex')) 

Jo Jx 



IX. APPENDIX C 



Here we provide exact bounds on 1\ and I2 given by (37 38). 
In order to estimate I2 we introduce the linear operator A : 



2x J x 2x\J~x Jq 

defined on the Hilbert space of the functions on [0, 1] with the scalar product: 



(f\g) = / x^x f(x)g(x)da 
Jo 



Then A is a bounded symmetric operator and: 



1 f (ab\Aab) 



(60) 



> 4e 2 \/2^[lne- 1 +0(1)] / dx I dx'g(ex,ex')y/x (61) 

JO Jx 

And we have: 



(62) 



Af{x) = ±- [ f(y)dy + — L= I f(y)Vydy (63) 



(64) 



(65) 



The unitary operator f(x) — > g(?/) = /(e y )e 5 / 4y from £ 2 ([0, l],xy^dx) onto £ 2 ([0, +00], dx) maps the operator A 
onto the operator A: 

-x/A rx x /4 r+00 

Ag(x) = — jf ey'*g(y)dy+ — \ e^g^dy (66) 

Then 



J i/cx _ \ ikx r ~x/A fr? 7 \ 

4(1/16 + fc 2 ) l/2 + *2fc [ } 
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Thus setting: 



9k(x) 



1 



|l + i4fc| 

{9k}k>o is a full set of pseudo-eigenvectors satisfying: 



[(1 + Uky kx - (1 - iAk)e~ ikx ] 



(68) 



4(1/16 + k 2 ) 



9k 



Thus the spectrum of A is (0, 4) and the spectral measure is purely absolutely continuous; the largest spectral value is 
4 with a pseudo-eigenvector g±{x) = x-\-4 corresponding to f±{x) — x _5 / 4 (4 — lnx). But ||/4|| is infinite and f± diverges 
at 0. The next step is to choose a family of functions b v such that a v = y/1 — b v is defined and {a^b^Aa^b^) /||^|| 2 
is close to 4. 

Thus setting f v (x) = mm(/4(x), 74(77)) for < 77 <C 1, we have: 



WfvW 2 
(fM) 



In 3 77 — — In 2 77 + 0(ln 77) 



In 3 77 - — In 2 77 + 0(ln 77) 
5 



Then: 



(fnWn) 

\w 



20 
I In 77 1 



+ 0(ln" 2 77) 



Thus frj is a good candidate for the linear part of the problem. Now, by the simple scaling: 



brj(x) 



fri(x) 



(69) 
(70) 

(71) 
(72) 



we get the nonlinear candidate satisfying b v (x) < 1/V% = y 1 — b 2 is well defined, a v (x) > l/y/2 and b v satisfies 
®- 

We must now estimate the simultaneous convergence of 12/11 ( [65] ) and /2 as 77 decreases. 

(m^) - (brjarjlAbqCLrj) = - (b v - b^a^A^ - b v a v } +2(b v - b^a^A^) 

< 2 (b v - brjdrjlAlbr,) 

= S(b v - brjCLrjlbrj) + 2 (6^ - 6^1 (A - 4)6^) 

< 8 - 6^6^) + 2||^ - b v aJ\\(A - 4)6J 



< 8 (brj - brjCLrjlbrj) + 8^ - b^UUb^A - 4)6^ 



where: 



(brj - brja^brj) = / 
JO 



b r] (x) 2 [1 — arj(x)] x\fxdx 



< 



/ b v (x) 2 [l-a r 
Jo 



(x)] 2 x\/xdx sup 



1 - a ri 



y/2 

< \\b v - 6 r? a r7 || 2 



y/2- 1 



and since b^ — b^a^ > 



ll&rj - MrjH HIM - ll^sll -2(6^-6^6^) 

<IIM 2 HM 2 
= Kll 2 
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By a direct computation: 



KII 2 <IIM 2 



5 1 In 77 | 



for 77 small enough, and thus: 



(brjlAbrj) (brjdrjlAbrjdrj) 8 ( 6 



| In ?7 1 \ V2 ~ 1 5 



2V6 



And finally, from (71) for b v and (73) and I\ — 1 1 fo^ 1 1 , (65) gives: 



I2<^- V 5 / 2 \\n 3 V \ 



+ 0(ln -2 rj) 

i + o- 1 



In 77I 



where 



C = 20 + 8 



a/2 6 



72- 1 5 



2V6 « 92 



(73) 

(74) 
(75) 

(76) 



Choosing 77 sufficie ntly small, this proves that AE 1 is strictly negative for any r s > 0; furthermore choosing 77 to 
minimize AE (Eq. 41), i.e., 9r s | In 77 1 2 = 2C7T, we obtain as r s goes to 0: 



AE < —r a exp 



5tt_ 0(1) 
3r. s A /r 



(77) 
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